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Abstract 



Hadronic processes in which hard diffractive production takes place have 
been observed and analyzed in collider experiments for several years. The 
experimental rates of diffractive W's and dijets measured at the Tevatron 
and the cross sections of diffractively produced dijets recently obtained at 
the HERA experiment are the object of this analysis. We use the Pomeron 
structure function obtained from the HERA data by two different approaches 
to calculate the rates and cross sections for these processes. The comparison 
of theoretical predictions with the measured values reveals some discrepancies 
that make evident conceptual difficulties with such approaches. A new version 
of the Ingelman-Schlein model is proposed as an attempt to overcome such 
difficulties and make theory and data compatible. 



PACS: 12.40.Nn, 13.60.Hb, 13.60.-r, 13.85. Qk, 13.87.-a 
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I. INTRODUCTION 



The phenomenological analysis of hard diffractive processes has become one of the most 
interesting theoretical laboratories to investigate the nature and structure of the Pomeron. 
The concept of Pomeron structure function was introduced by Ingelman and Schlein [Q] as 
an ansatz to investigate the eventual production of high-p^ jets in diffractive hadron in- 
teractions. Such a theoretical speculation has become reality when the UA8 Collaboration 
obtained the first measurements of diffractively produced dijets 0. However, further quan- 
titative analyses carried out at the CERN pp Collider have shown that the predicted 
rates obtained by the Ingelman-Schlein (IS) model were too much high in comparison with 
the measured values. 

The same problem appeared when similar measurements were performed at the Fermilab 
Tevatron Collider. The ratio of diffractive to non-diffractive dijets predicted by the IS model 
1^ for the Tevatron energy resulted to be more than one order of magnitude above the actual 
measurements. 

Despite the fact that the IS predictions about the existence of diffractive jets turned out 
to be correct, one might think that this model is acceptable only in a "qualitative" sense 
since the predicted rates were completely off the values actually measured. Let us examine 
this issue in a more detailed way. 

The calculational scheme proposed by such a model makes use of the so-called factoriza- 
tion hypothesis, that is in a reaction like p + p — > p+jets + X , the whole process is supposed 
to occur as a sequence of two independent steps: first a proton emits a Pomeron, then partons 
of the Pomeron interact with partons of the other proton producing (for instance) jets. This 
picture seems to be a sort of straightforward extension of the parton model to diffractive 
processes. In fact, it is quite appealing as such, but it is not obvious that the factorization 
property should apply to this type of hadron interactions (see and references therein). 

A similar kind of problem affects also another class of processes. Experimental results 
apparently in favor of factorization have come up from the DESY ep Collider, the so-called 
HERA experiment. Events of deep inelastic scattering (DIS) tagged with rapidity gaps, 
observed and analyzed by the HI and ZEUS collaborations [^,0 , have shown the same char- 
acteristic pattern exhibited by hadron diffractive dissociation processes in the kinematical 
region where Pomeron exchange is dominant. These observations strongly suggested that, 
in fact, the internal structure of the Pomeron was being probed. Moreover, this pattern, 
which resembles the peak observed in soft diffractive dissociation, seemed to be independent 
on kinematical variables other than x-p (the fraction of the proton momentum carried by 
the Pomeron), allowing that estimations of the Pomeron intercept were done 

These were the first results of diffractive structure function. More recently, further 
and more precise measurements performed over an extended kinematical region have shown 
evidences of factorization breaking in the xsp variable Such an effect, however, can 
be understood by an adequate computation of secondary- reggeon contributions 0, leaving 
open the possibility that factorization applies when Pomeron exchange is the dominant 
mechanism. 

Nevertheless, this conclusion has to be taken with some caution because other forms 
of factorization breaking, which are not evident from the xjp dependence, can occur. For 
the moment, let us assume (as the IS model does) that the factorization hypothesis is an 
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acceptable ansatz. This assumption leads us to more definite theoretical questions which 
were present in the IS model from the very beginning. 

From a quantitative point-of-view, factorization appears in the IS model as the product of 
two quantities representing the two-step process mentioned above: (1) the so-called Pomeron 
flux factor, which is supposed to give the probability of a Pomeron being emitted by a proton 
(or antiproton), and (2) the elementary cross section resulting from the interaction among 
partons belonging to the Pomeron and to the other proton [|I|. In order to calculate this 
elementary cross section, the knowlegde of the Pomeron structure function is, of course, an 
indispensable requisite. By the time the IS model was proposed there was no experimental 
information about the Pomeron structure and, thus, only estimations based on "educated 
guesses" could be done fl],^. Today, the HERA data of diffractive DIS are available and 
one can try to extract from them the Pomeron structure function. However, this is not a 
straightforward procedure because the results are dependent on the model one chooses for 
the Pomeron flux factor. 

In a recent paper 0, we have presented a study on the Pomeron structure function in 
which two different forms of flux factor were employed, one derived from the standard Regge 
theory |10| and the other obtained from the so-called renormalization procedure . The 



former assumes factorization whereas the latter implies in a sort of factorization breaking. 
We have shown that the quark/gluon content of the Pomeron changes significantly whether 
one chooses the standard or the renormalized flux factor. In the present paper, we apply 
these results to estimate the rates of diffractively produced ly s and jets and compare such 
estimations with the available data. 

As we shall see, no one of these models, standard or renormalized, is able to offer a 
completely satisfactory description of the data, although the latter is partially successfull. 
In order to overcome these difficulties, we introduce a new version of the IS model, which is 
quite intuitive and presents promising results. 

Although the comparison to data is central to our analysis, we emphasize that no attempt 
to fit the theoretical outcomes to the experimental rates or cross sections has been donea. 
In fact, such an attempt could conceal the problems that we intend to make evident. 

The paper is organized as follows. In Sec. II, we summarize our procedure to determine 
the Pomeron structure function from HERA data and present the parametrizations obtained 
in P that are used in this paper. In Sec. HI, we present the formalism used to calculate the 
cross sections for diffractive production of W^'s and jets. In Sec. IV, the results obtained with 
the standard and renormalized flux factors are shown in comparison with the experimental 
data. In Sec. V, we present a discussion about the IS model and a new approach is suggested. 
Our concluding remarks are in Sec. VI. 

II. POMERON STRUCTURE FUNCTION FROM HERA DATA 

The cross section for diffractive DIS processes. 



^We refer the reader interested in appreciating an analysis that tries to conciliate the IS model 
with data by fitting, to the paper of Ref . p . 
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e{k)+p{P) ^ e\k')+p'iP')+X{M'^), (1) 
is given by the expression 



dx dQ^ dxw dt x [ ^ 2[1 + x^, t)] 

where x, and y are the usual DIS variables. Besides these, two other variables are used 
to specify process ([l|), 

and 

3 = ^ ^ ^ (4) 

where xjp is the momentum fraction carried by the Pomeron emitted by the scattered proton 
and f3 is the fraction of the Pomeron momentum carried by the struck quark. In Eqs. (^) and 
(I), W is the energy in the 7*p center-of-mass system, t = {P — P'Y is the four-momentum 
transfer at the proton vertex and Mx is the invariant mass of the hadronic system X. From 
these equations, one obtains the relationship 

= ^- (5) 

The HERA data used to perform our analysis were obtained under the assump- 
tion that R^{x,Q^,XTp,t) = and, since t was not measured, the cross section must be 
considered as integrated over this variable, that is 

SO that the experimental data are expressed in terms of the dijfractive structure function 
F,''^'\p,Q^,xjp). 

These data have shown for the first time a very clear diffractive pattern, that is the 
characteristic xp-dependence observed in diffractive dissociation of hadrons. This feature 
was observed irrespective of the Q^) values considered and this was very suggestive that 
a factorized expression like 

F,''^'\p,Q^xjp)=g{xjp) F^{P,Q') (7) 

could apply!. Based on this factorization hypothesis and on the IS model it is usual to 
interpret g{x-sp) as the integrated-over-t Pomeron flux factor and F^{f3, Q"^) as the Pomeron 
structure function. 



^As mentioned in the Introduction, new data obtained by the HI Collaboration Q in more 
extensive kinematical regions of both (3 and have shown very clear factorization breaking which, 
however, can be understood as being an effect of secondary reggeon contributions other than the 
Pomeron. Our analysis |^ refers to data obtained in a particular kinematical region where Pomeron 
exchanges are supposed to be dominant. 
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Besides our analysis, there are others in the hterature that are based on a similar proce- 
dure (see for instance |0,12,1^). Our main concern, however, was confronting two different 
approaches: one in which the standard flux factor is employed and the other which corre- 
sponds to the renormalized flux factor (for brevity, we will refer to these quantities hereafter 
as STD and REN flux factors, respectively). For the former, it was assumed the Donnachie- 
Landshoff expression [ p!0|] , 

/s™(x:P,t) = ^ [Fi(t)]^x]p^"(^ (8) 



whereas the latter is determined from the procedure prescribed in [0, that is 

WxP,t)= (9) 

where 



N{xw^J = / ™^ dx^ / fsTD{xjP,t) dt. (10) 

J Xjp Jt= — oo 

By introducing Eq. (^) into Eq. (|1^) and assuming an exponential approximation for the 
form factor, -Ff (t) — e^°*^*\ one obtains 

N{x^^J =K^— [E,,(7 - 2e lna;np_) - E,,{^ - 2e Inxjp^J] (11) 

where Ei{x) is the exponential integral, K = dp^/in'^ and 7 = h^e/a' . 

An important point to be noted here is that, in the calculation of the soft diffractive 
dissociation cross section, the minimum value of x-^p is xjp^^^ = {nip + m^)^/s, but it is 
assumed that when one applies Eqs. (^-(|TT]) to the diffractive DIS analysis this quantity 
should become (see [0) 

XTP . (12) 

This distinction is pretty important and will be a matter of discussion in Section IV. 

With these flux factors (integrated over t) introduced into Eq. (^, the Pomeron structure 
function was obtained from F2^'-^''(/5, Q^, xp) data under the assumption that 

Ff{P,Q^)= E PW,Q^) + q0,Q^)] = '^pn^,Q^), (13) 

i=u,d,s 

with Il(/3, Q"^) = J2i=u,d,s [QiiP, Q^) + ^i(/5) Q^)] representing a quark singlet that would evolve 
with Q"^ according to the DGLAP equations fl^. The starting scale for the evolution was 
assumed to be Qq = 4 GeV^ and no sum rule was imposed on the parametrizations to 
perform the fitting (however, in one of the cases presented below, we use the sum rule to 
determine a parameter of normalization). 

Details about this procedure and about the results can be found in 0, but roughly 
speaking this analysis has put in evidence three major points: 



5 



• The quark/gluon content of Pomeron as obtained from HERA data via IS model 
depends strongly on which kind of flux factor is assumed. That means that the issue 
about (re)normalization is crucial. 

• STD flux factor favors a predominance of gluons at the initial scale of evolution 
while basically the contrary happens with the renormalized scheme. 

• Several different trials were carried out during the fitting procedure. In almost all of 
them, the initial quark and gluon distributions preferred a hard or super-hard shape. 
Thus, this analysis has practically ruled out the possibility for soft distributions at 
initial scale in diffractive DIS. 

For the present analysis, we have chosen two parametrizations for each flux factor. These 
parametrizations, the most representative of our analysis, are described below. 
Fit 1: These parametrizations were obtained in with STD flux and correspond to a 
combination that we call hard-hard, i.e. both quark and gluon distributions have a hard 
shape at the initial scale of evolution: 

/3S(/3,Q2) = 2.55/5 (1-/3), 

P giP,Ql) = 12.08 P (l-P). (14) 

Fit 2: This case refers also to parametrizations obtained with the STD flux; the initial 
distributions correspond to a super-hard profile imposed to gluons by a delta function while 
quarks were left free to change according to the data (we refer to this case as free-delta): 

/3S(/3,Q2) = i.5i po.6i (i_;5)0.84^ 

P g{P,Ql) = 2.06 6(1 -P). (15) 

Fit 3: These are parametrizations obtained with the REN flux factor and a initial combi- 
nation of the type hard-hard: 

p j:{P,qI) = 5.02 p (i-p), 

Pg{P,Ql) = 0.98 Pil-P). (16) 

Fit 4: These are also parametrizations obtained with the REN flux factor and a combination 
of the type free-zero, that is the quark distribution was left free while the initial gluon 
distribution was imposed to be null: 

/3S(/3,Q2) = 2.80 P'-^' (l-/3)°■'^ 

Pg{P,Ql) = 0. (17) 

In the case of Fit 3, it was very difficult to fix the normalization parameter for the gluon 
distribution (see [^]). In the expression used here, Eq. ([T^), this parameter was established 
by using the normalization obtained for quarks and imposing sum rule. In fact, it was 
because of such a difficulty that Fit 4 was performed. 

All these four combinations of flux factors and Pomeron structure functions were applied 
in the calculation that follows. 
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III. CROSS SECTIONS FOR DIFFRACTIVE PRODUCTION PROCESSES 



As we have seen in the previous section, several possible parametrizations for the Pomeron 
structure function are allowed by the HERA data depending on what one assumes for the flux 
factor. Our aim is applying these parametrizations to calculate the diffractive production 
rates of different processes in order to compare the results with the available data and 
analyze the implications. In this section, we present the cross sections used to perform these 
calculations. 

The generic cross section of a process in which partons of two hadrons, A and B, interact 
to produce jets (or VT's), 

A + B-^ Jets iW) + X, (18) 
is given by the standard parton model as 

d(TAB^Jets(W) = J2 fa/AiXa, fJ-"^) /fe/B (^^b, /i^) '^^"^-"f^(^) (19) 

a,b,c,d dt 

In order to adapt such a cross section to a hard diffractive interaction, one assumes 
(in the spirit of the IS model) that one of the hadrons, say hadron A, emits a Pomeron 
which is made up of partons itself. The procedure we adopt in such a case is replacing 
Xafa/A^Xa, fi'^) in Eq. ( |19|) by the convolution between the distribution of partons in the 
Pomeron, /5/a/]p(/5, /i^), and the "emission rate" of Pomerons by A, fw{xTp,t). The first 
quantity corresponds to distributions like those discussed in the previous section, that take 
part in the definition of the Pomeron structure function, while the second corresponds to 
the Pomeron flux factor. In such a case, we have 

Xa fa/A{Xa,fi^) = J dxjp J d(3 J dt /]p(X]p,t) (3 /a/]P (/5, /i^) 5{(3 - Xa/xw)- (20) 

Defining g{x-!p) = J^^dt fw{xw,t), one obtains 

Xa fa/A{Xa, fJ''^) = [ gixjp) Xa fa/wiXa/ X^ , fJ^"^) ■ (2l) 

J Xip 



This is the characteristic procedure that is applied below to calculate the cross section of 
diffractive processes. 



A. Diffractive Hadroproduction of 

In this analysis, we consider diffractive production in the reaction 

p + p^p+ e u)+ X, (22) 



which was experimentally studied at the Tevatron Collider by the CDF Collaboration [ITS 



We assume that a Pomeron emitted by a proton in the positive 2;-direction interacts with 
a p producing that subsequently decay into u. In this configuration, the detected 
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lepton (e"*" or e~) would appear boosted towards negative t] (rapidity) in coincidence with a 
rapidity gap in the right hemisphere. 

The cross section for the inclusive lepton production by this process isi 



= Ylj dxw g{xw) j dEr fa/w{xa, jJ''^) fb/p{xb, jJ''^ 

where 



GsTwMw 



Mw e^^ 



(23) 



(24) 



Xh 



Mw e~ 



and 



t = -Et Ml 



w 



A + JiA^-1] 



(25) 



(26) 



with A = Mw/'^Et- The upper signs in Eqs. ( |24D and (^) refer to production (that is, 
e"*" detection). The corresponding cross section for W~ is obtained by using the lower signs 
and t u (see the Appendix). 

Since the experimental data of diffractive W production presently available are not highly 
precise (in fact, CDF data [15] are the first measurements of this process), our calcula- 
tions consider only leading order contributions. In this sense, we note that studies of non- 
diffractive W production performed with CDF cuts indicate that next-to-leading (NLO) 
order corrections could be around 10% of the production cross section 



see 



T6||). However, 



the experimental information is given in terms of a ratio between diffractive to non-diffractive 
W production, thus the effect of NLO corrections in the final results are expected to be even 
smaller than this percentage. 

The contribution of other competitive processes like inclusive hadroproduction of W+jet 
and + 7 (and respective NLO corrections) are expected to be very small as well (about the 
former process, see for instance analyses by Giele et al. in Ref. and references therein, 
and the CDF account [15|, while with respect to the latter we refer the reader to the work 
of U. Baur et al. in Ref. 



B. Diffractive Hadroproduction of Dijets 

In order to calculate the process by which two hadrons, say p and p, interact generating 
dijets. 



P + P^ jl + J2 + X, 



(27) 



See a more detailed discussion about this cross section in the Appendix. 
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one starts from the expression 

dapp^j^j^ = fa/p{,Xa,ld^) dxa fb/pixb, fi^) dxb ^'^"^r'"^ di, (28) 

a,b,c,d '^^ 

in which it is assumed that partons a and b belong to hadrons p and p in their initial states, 
and partons c and d will give rise (in leading order) to the dijet pair in the final state. The 
distributions /a/p(a^a, Z^^) and fb/p{xb, f^"^) are the structure functions evolved to an adequate 
scale /x^ and daab~*cd/ dt stands for the QCD matrix elements proper to this calculation. 

With a suitable changing of variables, namely dxa dxb dt —>■ 2Et Xa Xb dEx drj' drj, one 
gets the differential cross section in terms of the rapidity t] of one of the jets. 



da ^ [^Tr^ax ^ rVLax „ „ d(J ab-. 

1- = 1^ "^T / dr] Xa Ja/p{Xa,IJ' ) Xb fb/p{Xb, ) -J" 



where 



= ^(6-^ + 6-"'), a;, = ^(e'' + e^'), (30) 

with Et being the jets transversal energy. These expressions apply to the non-diffractive 
case. 

By using the convolution procedure described above, the cross section for diffractive 
hadroproduction of dijets becomes 



-r = J2 [ dE^ I dr]' [ dxw g{xw) I3fa/-sp{f3,ii^) Xbfb/p{xb, n"^) 



where P = Xa/xw with Xa and Xb given by (|30|). The kinematical hmits for the above 
expression are 

In <ri'<\n (32) 

and 

E = ^ (33) 
with Et^.^, x-s>^.^, and xjp^^^ established by experimental cuts. 

C. Diffractive Photoproduction of Dijets 

The process considered now is diffractive photoproduction of dijets, obtained from the 
reaction 

e++j9 ^ e++p' + X(ji+j2 + X'), (34) 
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in which the positron is scattered at very small angles, implying that the emitted photon 
has a very low momentum (Q^ ~ 0) and can be considered, in a good approximation, as 
real. In such a context, the positron acts just as a source of photons that are emitted with 
a certain energy spectrum. This emission can be expressed in terms of a "photon flux" 
through the so-called Equivalent Photon Approximation (EPA) (or by the best known 
Weizsacker- Williams approximation) . 

There are experimental evidences that photoproduction processes take place by two 
mechanisms whose calculation is considered in terms of: (1) the direct component, in which 
the coupling of the photon to partons of the proton is point-like; and (2) the resolved com- 
ponent, in which the photon fluctuates to a partonic structure whose constituents interact 
with the partons of the proton. In fact, photoproduction of dijets is one of main processes 
by which the photon structure function is measured because this reaction is quite sensible 
to the quark/gluon content of the photon, even in leading order. 

In the case of diffractive photoproduction, according to the IS model, it is the partonic 
structure of the Pomeron that is probed by the photon itself (in the direct process) or that 
is envolved in interactions with photon constituents (in the resolved process). 

At this point, it is interesting to note that the way by which the IS model was 
conceived to describe hard diffractive production is in complete analogy to the EPA in 
photoproduction processes. Just as the electron (or positron) in photoproduction, the proton 
in a diffractive interaction is scattered at very small angles and practically does not take part 
in the effective reaction. In a analogous way to the emission of photons and to the idea of 
(equivalent) photon flux, it sounds natural to talk about Pomeron emission and "Pomeron 
flux factor". However, the problem with this analogy is that while the photon flux has a 
well-based theoretical derivation in QED, the Pomeron flux factor is obtained in a totally 
phenomenological way. This issue is central and will motivate more discussion below. 

Back to photoproduction, the momentum distribution of the interacting object that 
comes from the positron vertex (namely, the photon itself or its constituents) is given by 

Xafa/e{xa,H^) = J dQ^ J dxj J dy G{y,Q'^) x^fa/^{x^,i/) 5{x^ -Xa/y), (35) 

where the photon emission is described by the flux G{y,Q'^), obtained in the EPA context, 
and x^fa/^{x^, fi^) is the photon structure function, with x^ being the fraction of the photon 
momentum carried by partons. 

The derivation of the photon flux can be found elsewhere (see, for instance, 
integrated form reads 

Giy) ^ fdQ'G{y,Q') ^ \[l + [l - yf]\n^ - 2{l - y)] , (36) 

with Q^^j^, given by the experiment and Q^m — ^1 2/^/(1 ~ y)i so that Eq. (^) can be 
written as 

Xafa/e{Xa, fJ''^) = dy G{y) — fa/j{ — , jJ,'^) . (37) 

J y y 



19 ) and its 
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1. Cross section for the resolved component 



The cross section for diffractive photoproduction relative to the resolved component can 
be obtained in an analogous way to the hadroproduction expression, Eq. (0), but using 

xa = yx^ = ^{e^ + e^'), Xi, = (3x^ = ^{e'^ + e""'), (38) 

and Eq. (^) so that 

^ = y dE^ J dri' j dy G{y) j dxw g{xw) [3 fb/-sp{[3, j/) Xyfa/^{x^,i/) (39) 

with integration limits established likewise (the limits for the variable y are given by the 
experiment). 

Diffractive photoproduction data are also given in terms of other cross sections, da/dEx, 
da/dW, da /dp and da/dx^^ whose explicit expressions can be obtained from Eq. ( p9| ) with 
some appropriate change of variables. 



2. Cross section for the direct component 

The cross section corresponding to the direct component is obtained just by replacing 
the photon structure function by 

fa/^iXy,f/) = 6{1 - X-y) (40) 

in Eq. (|39|) and in the other cross sections for the resolved component. 

Note that, in this case, the Bjorken variable for the elementary process in leading order 
is Xa = y- Another consequence is that there is no direct component for da/dx-y in leading 
order. This will have important implications in the comparison of theoretical calculations 
with experimental data. 



IV. RESULTS AND DISCUSSION 

In the following, we present the results of our calculations of hard diffractive production 
processes whose cross sections were discussed in the previous section. For all of them, we 
have considered the four possibilities of Pomeron structure function discussed in Section II. 
As for the proton and photon structure functions, we have used the GRV (Leading Order) 
parametrizations pO| , ^ . 



A. Diffractive T^'s and the CDF rate 

As mentioned before, the diffractive W production has been studied at the Tevatron 
Collider {^/s = 1.8 TeV) by the CDF Collaboration in terms of a particular decay 



mode, that is W ^ e v. The measurements were triggered by electrons (or positrons) with 
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transversal energy Et > 20 GeV in the central region, \ri\ < 1.1, and corresponding to 
xjp < 0.1. 

In Fig. 1, we show the results obtained with Eqs. (p3D-(p^) for the STD flux factor (upper 
part) and for the REN flux factor (lower part). The calculations were performed with the 
CDF kinematical inputs, assuming that the Pomeron is emitted by a proton directed towards 
the right hand side. Due to this fact the distributions are boosted towards negative 
rapidity leaving an empty space (the characteristic rapidity gap) for rj > 1.5. In the same 
figure, we show the cross section for non-diffractive W production for comparison. 

It is clearly seen in Fig. 1 that the diffractive W production is more abundant for the 
STD flux than it is for the renormalized one in spite of the fact that the Pomeron structure 
function applied to the latter case is much richer in quarks (see Section II). 

In Table |, we present the ratios of diffractive to non-diffractive production rates calcu- 
lated with the different combinations of flux factor with IP structure function in comparison 
with the experimental value. The theoretical ratios were obtained by integrating the cross 
sections shown in Fig. 1 over the CDF limits, — 1.1<?7<1.1. As can be seen, the cal- 
culations with the STD flux factor overestimate the experimental rate by factor around 3. 
The results obtained with the REN flux factor, on the contrary, are much closer to the 
experimental value although a little below. 



B. Diffractive jets and Tevatron data 



Diffractive dijets rates were measured at the Tevatron Collider by the CDF p2| , p3| and the 
DO collaborations . In the CDF measurements two procedures have been used to isolate 
the diffractive events, one by the rapidity gap technique ||22[ and the other by using roman 



pots to detect the recoiling proton (or proton remnant) [Q. The DO analysis was performed 
only by using the rapidity gap technique, but it includes rates for two energies, a/s = 630 
GeV and y/s = 1800 GeV ||2^. In Table ||, the rates obtained by these experiments are 



shown as well as the kinematical cuts implemented in each case. We should notice, however, 
that the CDF rate given in the column (a) is the only experimental value already published, 
the others have to be taken as preliminary results. 

In Fig. 2 (a)-(d), we show the inclusive diffractive jet cross section calculated with the 
Pomeron structure fuctions given by Fits 1-4 (the kinematical cuts corresponding to each 
figure is identified by the letter in the top of the columns of Table II). The non-diffractive 
jet cross section calculated with the respective kinematical cuts is also shown. 

A characteristic feature of these calculations is that the results obtained with the STD 
flux with both fits are much higher than those given by the REN flux, whose fits in turn 
produce rates practically indistinguishable. 

From these curves, we have the ratios of diffractive to non-diffractive production rates 
given in Table |IT1| to be compared with the experimental values. Again the values obtained 
with STD flux are much larger than the actual measurements while the results given by the 
REN flux are in general agreement with the experiments. 
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C. Diffractive jets and ZEUS data 



The experimental data of diffractive photoproduction of jets used in this analysis were 



obtained by the ZEUS Collaboration at the HERA experiment [^], with the energy of 
the 7*p system between the limits 134 < W < 277 GeV and with the photon virtuality 
restricted by < 4 GeV^. Other kinematical variables that specify the outcomes of this 
experiment are the following: -1.5 < rj^jf^ < 1, E^^^^ > 6 GeV, and 0.001 < xsp < 0.03. 
Due to the asymmetry between the positron and proton beam energies, E^'^'^ = 27.5 GeV 
and E^p'' = 820 GeV (which correspond to a/s = 300 GeV), the rapidity variable in the 
center-of-mass system (cms) is given by r]^^^ = rf'^^ + (1/2) \n.{E'-^'^ / E^^^) . 

As mentioned before, the experimental differential cross sections are given in terms of 
da/drjjet, da/dW, da/dEx, da/dp°^^, and da/dx'^'^ . The superscript in the variables (3°^'^ and 
x°^^ indicates that these are quantities not directly measurable, but instead are observables 
obtained from the jets kinematics (see details in [^). 

In Fig. 3, we show the results of da/drjjet for both STD and REN flux factors. Here 
appears a situation partially different from what we have seen in the previous cases: the 
results obtained with the STD flux continues to overestimate the measured cross section by 
a large extent, but those obtained by the REN flux are not compatible with the data as it 
was previously observed in diffractive hadroproduction of W^s and jets. 

Basically the same features are seen in Figs. 4(a)-4(d), in which we present da/dW, 
da/dEx, da/djS"^'^, and da/dx°^'^. It is important to notice, however, that the curve whose 
shape most correspond to the /^-distribution in Fig. 4(c) is that one obtained with super- 
hard gluons (Fit 2). An additional observation about the results of Fig. 4(d) is that the 
cross section for this case does not include the direct component since our calculation are 
performed only to leading order. It is known, however, that in next-to-leading order the 



direct component presents appreciable contribution for x^ > 0.75 pBj 



D. Discussion 



Taking as whole, the results presented in Tables I and HI, and in Figs. 1-4 show that the 
combination STD flux factor with the Pomeron structure functions obtained from Fits 1 and 
2 cannot describe the diffractive production of I^'s and jets discussed here. These results 
are systematically above the data, sometimes by an order of magnitude or more. Thus, a 
natural conclusion seems to be that the IS model with the STD flux factor is ruled out by 
the experimental data. 

The other general observation is that the results obtained with the REN flux factor and 
Fits 3 and 4 are in a reasonable agreement with the diffractive hadroproduction data, but 
the same theoretical scheme fails to describe diffractive photoproduction. 

Since the renormalization procedure is usually presented as a way of reconciling the 
Ingelman-Schlein approach with the experimental observation (and, in fact, it seems to work 
in hadroproduction), the question is why such failure happens in diffractive photoproduction. 
The explanation is in the renormalization factor itself and is given in the following. 

As pointed out in Section II, for the case of diffractive DIS, the lower limit that enters 
in the definition of the (re)normalization term, Eq. (p!OD, is xjp^-^^ = Q'^/j3s. In diffractive 



13 



photoproduction, we have a (dynamically) similar process, except for the range of values 
assumed by Q^, implying that the definition of xp,^.^ should be the same. In order to be 
clearer about the way the Pomeron structure function was obtained in the renormalized 
case, we rewrite Eq. as 

F,''^'\p,Q',xjp) = gREN{xw,xjp^J F^iP,Q'), (41) 

meaning that, in such a case, the Q^-dependence comes from the DGLAP evolution and 
from Xsp^-^ implicit in the renormalization factor. 

In all of the cross section calculations presented above (including, of course, the diffrac- 
tive ones), we have applied the procedure usual in QCD/parton model of using Et as the 
evolution scale in the structure functions. However, in order to be consistent with the 
original proposal |jll|, in diffractive photoproduction calculations we must assign to the Q^- 



dependence that belongs to renormalization factor, the Q -values referring to the ZEUS 



experiment 25 



In fact, there is a problem about which value to choose since < 4 GeV^ with the 
median value in ^ 10"^ GeV^ |2^. Even in disagreement with the data, the curves 



shown in Figs. 3 and 4 represent the most favorable situation, which corresponds to put 
= 4 GeV^. If one applies the median Q^, these results are reduced to a really small 
fraction, five times lower than the curves presented. This is easily understood from the fact 
that ~ in this kind of experiment and the smaller is, the larger the normalization 
factor becomes, reducing considerably the calculated cross section. 



V. HARD DIFFRACTION: CLUES TO A NEW APPROACH 
A. Diffractive Parton Model 

We propose here a new version of the Ingelman-Schlein model that, in our view, seems 
to be able of overcoming the difficulties presented and discussed above. 

First of all, we would like to state that the Pomeron flux factor, as it is presently es- 
tablished, is an ill-defined and misleading quantity that cannot be supported only by the 
analogy with the photon fiux factor (which seems to be best justification for it) or something 
alike. The concept is interesting, but its definition in terms of the (standard) Triple Pomeron 
Model leads to wrong results (as we have shown). Maybe, in the future, QCD will provide 
a rigorous definition for the Pomeron fiux factor but, at the moment, we see no reason to 
keep it. 

We think that, starting from the idea that the Pomeron is constituted of quarks and 
gluons, what one really needs to estimate the measured hard diffraction cross sections is 
a probability distribution that would connect hard interactions that occur at partonic level 
to the hadronic level at which diffractive processes are detected. We propose that such a 
distribution be given by the normalized function 

a^/ dt dxsp 
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where a^^/ represents the single diffractive cross section integrated over only one hemisphere. 
The other term is, of course, the differential cross section, that we assume to be known and 
that is, in principle, in agreement with the experimental data (that is what the superscript 
exp means). 

Let us call this quantity defined by Eq. p^), FsAxw-,'^-, diffraction factor since it rep- 
resents the probability distribution that a diffractive interaction takes place. Once the 
diffraction factor is known, we propose that the cross section of hard diffraction processes is 
the result of the convolution product 



dx-^ dt 



I = Fsd{xTp,t) (g) Tjhard, (43) 
HD 



in which T^hard stands for all the elementary cross sections involved in the specific process 
under consideration. Operationally, Eq. (|43|) represents what has already been done in 
Section III, if one replaces the Pomeron flux factor by the diffraction factor here introduced, 
since the convolution product is conceived to be taken in the same sense of Eq. (|20|). 

Eq. (1^ ) is, of course, reminiscent of the IS expression, with exception of the normaliza- 
tion that, in this case, is established by a'^Y instead of the Pomeron-proton cross section, 
ajpp, which is the original assumption (Cf. |I|). However, this small change implies in two 
important differences: 1) -Fsrf(x]p, t) is a normalized distribution by construction, and 2) 
is an experimentally observable quantity while ajpp is a model dependent one. In order to 
have a brief form to refer to it, we are going to call the combination of Eq. (^) with Eq. (^2]) 
Diffractive Parton Model (DIFFPM). 

Now, we intend to show that the renormalized flux factor is nothing but an approximate 
expression for the diffraction factor Fgdixjp, t). In order to do that, let us turn our attention 
to the single diffractive cross section as it is given by the standard Regge theory, 

^fi^ = fsMx^,t) a^p, (44) 

which is known from long ago to violate unitarity. Let us assume that the corrections 
necessary to make this cross section compatible with the data are (in principle) known from 
some physical effects (screening corrections, flavoring, etc.) and that they can be represented 
by a function C{s), which depends only on the energy, such that 

dt dxjp C{s) dt dxjp 

Function C{s) must not be confused with the renormalization factor N{x-sp^.^) at this point. 

Implicit in Eq. (|45|) , there is the assumption that the xjp and t dependences given by the 
STD model are in agreement with the data and that the real problem has to do only with 
the energy dependence. This assumption seems to be supported by the analysis presented 



m 



For simplicity of reasoning, let us momentarily assume that the Pomeron-proton cross 
section is constant, a^pp = ctq (which, in fact, it approximately is). From these hypotheses, 
we can extract two results: 
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Result 1: By replacing Eq. with (|4^) into Eq. one obtains 

Fsd{xw,t) = -x^ "0 — -, 46) 

/xff It=-oofsTD{xTP,t) dXjp dt 

that is the same expression of the renormalized flux factor, Eq. but in which it is imposed 
that xjp^^^ = {nip + rriT^y/s always, that is by deflnition. 
Result 2: Now, by integrating Eq. (^5]), one gets 

'^sT = 7^ / / fsTD{xjp,t) dxw dt, (47) 



Cis) 



w 



rmn 



from which we see that C{s) will be the same as N{x-sp^-^) if (and only if) a^^^ = do, and 
that means a^Y constant. From this reasoning, one can obtain the renormalized expression 
for soft diffraction in two steps: flrstly, one assumes that a^^^ = ctq and determines C(s) 
from Eq. (|17|), and secondly, one replaces C(s) so obtained in Eq. (^). In the resultant 
expression, uipp is not assumed to be constant anymore, but a constant factor, that is part 
of it, is adjusted according to the ISR data (Cf. JTTH). 

Let us discuss these results, starting from the second. By looking at the energy depen- 
dence of the datai, we see two different behaviors: from low energies up to the ISR 
energies, the cross section is clearly increasing, but from the ISR to the Tevatron energies, 
it is pratically constant (although data are really scarce in such a region). The latter is the 
region in which the renormalization scheme is applied, which is consistent with the above 
argumentation. In fact, the energy dependence obtained for the renormalized cross section 
is really mild, changing very little over a range of practically 10^ GeV (see Fig. 1 of p^). 

Of course, all of this is valid only under the supposition that follows this almost 
constant trend also in the empty region between the ISR and the Tevatron data. If that 
is not true and the cross section have some strong variations in this non-observed region 
and/or beyond the Tevatron energy, then the renormalization scheme is not valid for soft 
diffraction anymore. In this case, it would be necessary a function C{s) different from the 
renormalization factor iV(a;ip^-„), that would represent such variations. But, independently 
on that function, the diffraction factor given by Eq. (^6|) (in other words, the renormalized 
flux factor) would remain the same. Therefore, so far the conclusion is: even if the renormal- 
ization procedure were not the correct solution for the unitarization of the soft diffractive 
cross section, the renormalized flux factor would remain valid as an approximate expression 
for Fsd{x-!p,t). 

B. Application to diffractive photoproduction 

The preceding discussion allows us to change the line of argumentation and the way of 
looking at the theoretical results presented here so as to put in evidence the Diffractive Par- 
ton Model given by Eqs. (H2|) and (|13|). From this point-of-view, we have already shown that 



^We remind the reader that the experimental data of single diffractive cross section are conven- 
tionally established as (Tg^j = 2 crlY- 
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the DIFFPM was able to give a reasonable description of the diffractive hadroproduction of 
W^s and jets through an approximate expression for the diffraction factor given by Eq. (^6|), 
or in other words, by the renormalized flux factor. 

Now we are going to show that, despite the difficulties pointed out previously, it is possible 
to give a reasonable description for diffractive photopro duct ion with the same parametriza- 
tions for the Pomeron structure function obtained with REN flux factor in 0, but by 
applying the DIFFPM. In order to explain how that is possible, we need to consider the 
renormalization procedure again. For the sake of simplicity, instead of using the full ex- 
pression given by Eq. (0), let us take for iV(a;ip^-^) the approximate formulas given in 



Ref. Ill 



N'^'ix^^J ^ c, s'^ (48) 

and 



N^'^ix^^J^c^ (^1 , (49) 




which refer to soft and hard diffraction, respectively (ci and C2 are just constant factors). 
Based on the above equations, we see that an approximate form for Eq. ( ^T]) is 

F°<''(A(3^.p)«-|g^;^^^•-(ft(^^), (so) 



which can be rewritten as 



gSTD[Xw) 



'^)'Vj^(/3,g^) 



(51) 



except for constant factors. 

Therefore, we see that the term above between square brackets can be reinterpreted as an 
effective parametrization for the Pomeron structure function, in which the Q^-dependence 
comes from both the renormalization factor and the DGLAP equations. The term on the 
left, gsTo/N^^ , represents nothing but the diffraction factor as it is given by Eq. (^61) 
(integrated over t, of course). 

The above reasoning is a simplification to understand what has actually been done. In 
summary, we have considered the Q^-dependence that comes from the renormalization factor 
as part of the Pomeron structure function and, as such, it has worked as the evolution scale 
in the photoproduction calculations as well. By doing so, we have effectively established 
the unique renormalization factor, which is consistent with the DIFFPM. 

By applying this procedure to the formalism described in Section III-C, we have obtained 
the results for diffractive photoproduction of dijets shown in Figs. 5 and 6. In Fig. 5, we show 
the results for da/drjjet in comparison with the ZEUS data. As explained in Section III, these 
results are obtained by summing up the direct and the resolved components. To illustrate 
the importance of taking into account both contributions, we show these components for 
the dashed curve as a hachured area (direct component) and as a shaded area (resolved 



17 



component). Thus, now the theoretical results are compatible with the data and there is no 
ambiguity whatsoever about how to treat the Q^-dependence. 

In Figs. 6(a)-6(d), we show again the results for da/dW, da/dEx, da/dp°^'^, and da/dx°^'', 
but now obtained with this new procedure of calculation. From Figs. 6(a) and 6(b), we 
could say again that compatibility with data is achieved (less for 6(b)), but the same does 
not happen for the data of Figs. 6(c) and 6(d). In order to describe the data exhibited in 
Fig. 6(c), a super-hard distribution for gluons seems to be indispensable. Such a distribution 
is very likely to affect also the results shown in Fig. 6(d), that would tend to become harder, 
providing a better description for the data. However, as mentioned earlier, next-to-leading 
order calculation for the direct component would be necessary. 

Even not obtainning a perfect description of the data, we think that the combination 
of these results (shown in Figs. 5 and 6) with the results previously shown for diffractive 
hadroproduction composes a picture from which it is possible to appreciate the possibilities 
of the DIFFPM. 



VI. CONCLUDING REMARKS 

We have presented in this paper an analysis on hard diffractive processes with the aim 
of describing the W and dijet diffractive production rates measured at the HERA and 
Tevatron colliders. In order to do that, we suggested a new approach which represents a 
sort of modification in the Ingelman-Schlein model. 

One (obvious) weak point in the discussion presented in Section V-A is that we do not 
make any attempt to determine the function C{s) without the renormalization scheme. Fur- 
thermore, we assumed that all corrections could be concentrated in this factorized function 
which would depend only on the energy. Of course, this represents an over-simplification of 
what really happens, but the spirit was to show by a sort of toy-model that, even not know- 
ing everything about the physics of these corrections, it is possible to find some acceptable 
justification for the renormalized flux factor within the DIFFPM scheme. 

The real solution might be (it certainly is) something much more elaborate like, for 
instance. Tan's "flavoring" model [^, the "damping factor" proposed by Erhan and Schlein 
||29|| , or the "screening corrections" of Gotsman, Levin and Maor |^0| (or even something 
else). Whatever is the "right" solution for the problem of soft diffraction unitarization, the 
calculational scheme represented by Eqs. (|4^) and (|43|) would remain basically the same 
since it does not depend on the particular model used to describe the single diffractive cross 
section (as long as such a model be able to provide a good description of the experimental 
data). 

The basic idea underlying our proposal is that the probability of having a diffractive (= 
rapidity gap) event in soft or hard diffraction can be represented, in a good approximation, 
by the same function. In our approach, it is given by the diffraction factor here defined. 
We have shown that, with this idea, it is possible to give an acceptable description of 
the experimental data for both hadroproduction and photoproduction in hard diffractive 
interactions. 



18 



ACKNOWLEDGMENTS 



We would like to thank the Brazilian governmental agencies CNPq and FAPESP for 
their financial support. 

APPENDIX A: DIFFRACTIVE HADROPRODUCTION OF 
1. Non-difFractive W production 



In order to calculate the cross section for the reaction (|22|), we start by dealing with the 
expression for the non-diffractive production of W's, 

da 
di 



da = J2[dXaf dxb fa/p{xa) fb/p{xb) -^io-b ^ eu) di. (Al) 

a,b J J dt 



In this equation, fi/A{xi) refers to the distribution of partons i in hadron A. In the hadron 
center-of-mass system (cms), the total, longitudinal and transversal energy of the electron 
(or positron) are given, respectively, by 

Ee = ^ [Xa{l+COSe)+Xb{l-COSe)], (A2) 

= ^ [Xa{l + C0S9) - Xb{l - cos9)], (A3) 

Er = ^sn.e, (A4) 

in which the constraint s = XaX^s = has been used and where 6 is the electron scattering 
angle with respect to the proton beam direction (which, in this paper, is assumed to be the 
positive 2;-direction) . From the expression for the electron rapidity, 

V. = In (^^) . (A5) 



and using Eqs. ( |A2D -(|A^), one gets 



Vl + cos^ ^^^^ 



and 



XaS V S V 1 ~ COS 

The Mandelstam variables for the elementary process, a b ^ W ^ e are 

s = {pa+Pb)=M^, (A8) 

i=iPe-Pay = -^ (1-COS^), (A9) 
U=iPe-Pby = ~^ (l+COS^), (AlO) 
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and, therefore, we have dxh dt = xi^j^ — 1 dE^ drj^, where A is defined as 

Mw 



A 



2Er, 



Such a change of variables allows one to rewrite Eq. 



^ =J2 [ [ dXa fa/p{Xa) fb/piXb) 



- 1 dt 



as 

^(ab ^ W ^ eu) 



(All) 



(A12) 



Now, from Eqs. ([A4|) and Eq. ( |A11| ), one obtains 

cos 6 - 



Va^ 



A 



(A13) 



Of course, the positive or negative signs indicate the direction in which the electron (or 
positron) is being emitted. This sign is chosen according to the following criterion: due 
to helicity conservation, the electron is preferentially emitted in the proton beam direction, 
such that in the case of W~ — > e~ one should use 



cos 6' = 4 



A 

By applying such a criterion to the case of W~ production, from Eqs. 



(A14) 



one 



gets 



Xn 



Xb 



w 



^ ,A-^A^ - 1 



u 



[I + cos( 



-Et Mw{A + - 1 



-Et Mw {A - VA^^) . 



(A15) 

(A16) 

(A17) 
(A18) 



Similarly, in the reaction v^. the positron is preferentially produced in the 

antiproton direction, such that 



cos6' 



A 



Thus, in the case of production, Eqs. ([A^)-( |Ani| ) give 

(A^^fW^ 



Xa ~ 


_ Mw 


Xb = 


Mw 


u = 


-Et Mw 


i = 


—Et Mw 



A - - 1 



(A19) 

(A20) 

(A21) 

(A22) 
(A23) 
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The elementary cross section for W production is given by 



— - = S{Xa- XbS) Vab 

dv XbS bMwi- w 



(A24) 



where Mw is the W mass, Vab is the Kobayashi-Maskawa matrix element, is the Fermi 
constant, and Vw is the W decay width. In Eq. (|A24|) , variable v stands for m or t according 
to W~ or production, respectively. 

In summary, the cross section for W~ production is 



GsTwMw 



u 



(A25) 



for which Eqs. ( |A15| )- (|A18|) apply. For the case of production, the cross section is 



da 
drje 



J2 J dEr fa/p{Xa) fb/piXb) 
a,b 

with Eqs. (|A2UD-([A2^). 



(A26) 



2. DifFractive W production 



From the discussion above, the cross section for the diffractive case is easily obtained. 
Introducing the prescription established by Eqs. (|20| ) and ( pT] ) into the expressions derived 
above, the cross section for diffractive W production becomes 



da 
dr]e 



J2 J dxjp g{xjp) j dEr fa/w{xa) fb/p{xb) 

a,b 



GsTwMw 



(A27) 



where g{xTp) is the integrated flux factor and 

Mw 



Xa 



Xb 



X 

Mw e 



(A28) 
(A29) 



We remind that, in the definition of these variables, it is assumed that the Pomeron is 
emitted by the proton and that W production is the result of ff — p interaction (as in the 
CDF experiment). The choice of the signs and of the variables t and u proceeds in the same 
way as in the non-diffractive case. 
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FIGURES 




FIG. 1. Rapidity distributions of e emitted in W hadroproduction processes. The curves 
labeled ND refer to the non-diffractivc cross section. The other curves correspond to the result of 
diffractive production, with the labels indicating the respective Pomeron structure function and 
flux factor used in the calculation (see text) . The vertical dotted lines establish the rapidity limits 
within which the CDF measurements were performed. 
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FIG. 2. Rapidity distributions of jets from hadroproduction processes. The curves labeled ND 
refer to the non-diffractive cross section (multiplied by a factor 10~^). The other curves correspond 

to the result of diffractive production, with the labels indicating the respective Pomeron structure 
function and flux factor used in the calculation. The letter in the top of each figure indicates the 
corresponding kinematical cuts presented in Table II and applied to the calculations. 
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FIG. 3. Rapidity distributions of jets from diffractive photoproduction processes. The labels 
indicate the respective Pomeron structure function and flux factor used in the calculation. The 
experimental data were measured by the ZEUS Collaboration 
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FIG. 4. Cross sections relative to diffractive photoproduction of jets in comparison to ZEUS 
data 1 25]. As in the other figures, the curve labels indicate the respective Pomeron structure 
function and flux factor used in the calculation. 
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• data - ZEUS (1994) 
Fit 3 (DIFFPM) 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 



^Jet 

FIG. 5. Rapidity distributions of jets from diffractive photoproduction processes. The labels 
indicate the respective Pomeron structure function used in the calculation, but with the redefinition 
discussed in the text. The experimental data were measured by the ZEUS Collaboration For 
the dotted curve, it is also shown its components: the direct contribution (hachured area) and 
resolved contribution (shaded area). 
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TABLES 



EXPERIMENT 

(W's) 


RATE 


STD 
Fit 1 
(hard-hard) 


STD 
Fit 2 
(free-delta) 


REN 
Fit 3 
(hard-hard) 


REN 
Fit 4 
(free-zero) 


CDF (Rap-Gap) 


1.15 ±0.55 


3.12 


3.54 


0.53 


0.58 



TABLE I. Diffractive production rates of W's (all values are given in percentages). 



EXPS. 


(a) 
CDF 
(Rap-Gap) 


(b) 
CDF 
(Roman Pots) 


(c) 
DO 
(630 GeV) 


(d) 
DO 
(1800 GeV) 


RATES 

(%) 


0.75 ±0.10 


0.109 ±0.016 


1-2 


0.67 ±0.05 


rapidity 


-3.5 < r? < -1.8 


-3.5 <r]< -1.8 


-4.1 < r? < -1.6 


-4.1 < r? < -1.6 




xw < 0.1 


0.05 < Xff < 0.1 


xip < 0.1 


x-sp < 0.1 


Et • 


20 GeV 


10 GeV 


12 GeV 


12 GeV 



TABLE IL Experimental data and respective kinematical cuts for different measurements of diffrac- 
tive production of dijets. 
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EXPERIMENTS 
(JETS) 


RATES 


D LU 
Fit 1 
(hard-hard) 


D LU 

Fit 2 
(free-delta) 


Fit 3 
(hard-hard) 


Fit 4 
(free-zero) 


CDF (Rap-Gap) 


0.75 ± 0.10 


15.3 


6.33 


0.62 


0.52 


CDF (Roman Pots) 


0.109 ±0.016 


3.85 


1.13 


0.15 


0.16 


DO (630 GeV) 


1-2 


15.4 


6.41 


0.87 


0.71 


DO (1800 GeV) 


0.67 ±0.05 


16.6 


6.14 


0.65 


0.57 



TABLE III. Diffractive production rates of dijets (all values are given in percentages). 



30 



